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Measuring the current -voltage characteristic of organic bulk heteroj unction solar devices some- 
times reveals an s-shaped deformation. We qualitatively produce this behaviour by a numerical 
device simulation assuming a reduced surface recombination. Furthermore we show how to exper- 
imentally create these double diodes by applying an oxygen plasma etch on the indium tin oxide 
(ITO) anode. Restricted charge transport over material interfaces accumulates space charges and 
therefore creates s-shaped deformations. Finally we discuss the consequences of our findings for the 
open circuit voltage V oc . 
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I. INTRODUCTION 

Since the efforts to build organic photovoltaic devices 
began, groups around the world tried to create the ideal 
semiconductors and device configurations!^ The first ex- 
periments with bilayer solar cells where soon superseded 
by devices with a blended active layers of two different 
semiconductors. Those bulk-heterojunction solar cells 
nowadays achieve power conversion efficiencies of up to 
8 % bringing the economical break even within reachPEl 

Experiments with innovative organic semiconductors, 
but also well-known systems, sometimes show an s- 
shaped deformation of the device current -voltage charac- 
teristic. Instead of an exponential current-voltage ratio 
as expected for diode structures, the response of such 
a device in conducting direction shows a local satura- 
tion and a later again increasing current around at a 
certain applied voltage region. A similar behavior was 
reported for inorganic copper indium gallium diselenide 
(CIGS) based devices which was attributed to the influ- 
ence of a chargeable cadmium sulfide buffer layer, creat- 
ing an illumination dependent energetic barrierP In or- 
ganic multilayer solar cells this deformation was observed 
several times in literature and was attributed to surface 
dipolesPl 

For numeric simulations a careful description of device 
interfaces therefore should lead to a double diode charac- 
teristic. In most publications metal-organic interfaces are 
described by a Schottky contact. The boundary condi- 
tion for the charge carrier densities therefore is defined by 
the thermionic emission theory resulting in a fixed charge 
carrier density at the interfaceJS"^ Restrictions in the 
ability to transfer charges from one side of the junction 
to the other can be expressed by the surface recombina- 
tion ratepSl By the latter, models such as mirror charge 
effects at metal-organic interfaces can be described, re- 
sulting in surface charge carrier densities different than 
predicted by the thermionic emission theory. 13 For or- 
ganic solar cells the importance of low surface recombi- 
nation rates on the power conversion performance was 



already discussed, mainly for minority charges! 14 * 15 ! 

Beside these boundary conditions at the interfaces var- 
ious models for physical effects inside the solar cell bulks 
are available for numerical simulations. Optical interfer- 
ence calculations are used to determine the spatial charge 
carrier gen eratio n by absorbed light for different device 
geometries! 16 * 17 ^ Assuming the existence of charge trans- 
fer states during the generation process an electric field 
dependence can be assigned based on the assumption of 
Gaussian distributed exciton binding distances.^ From 
the energetic and spatial disorder in organic materials 
charge transport models for the mobility of electrons and 
holes could be deduced using Monte Carlo simulation 
techniques E2HH] Combinations of those models ar e of- 
ten used to compare and evaluate experimental dataP3221 
Nevertheless none of these models is able to predict s- 
shaped current-voltage characteristics in simple device 
geometries as discussed in this paper. 

Within this article, we will show how to experimentally 
create s-shaped current -voltage characteristics of organic 
bulk heteroj unction solar cells by an extended oxygen 
plasma etching of the indium tin oxide (ITO) anode. Us- 
ing a numerical simulation, we qualitatively reproduce 
this double diode assuming reduced surface recombina- 
tion velocities at the anode. We will discuss the effects 
able to create s-shaped deformations as well as the conse- 
quences of our findings for organic photovoltaic devices. 



II. METHODS 
A. Experimental methods 

For the solar cell preparation structured ITO glass 
was successively cleaned in soap water, aceton and iso- 
propanol in each case for 10 min in an ultrasonic bath. 
The anode was exposed to an oxygen plasma etching with 
a certain nitrogen concentration for various time scales 
to change the ITO work function.^ On this substrate, 
an optional layer of poly (3,4-ethylenedioxythiophene) 



poly(styrenesulfonate) (PEDOT:PSS) was spincoated 
and transfered into a nitrogen atmosphere before anneal- 
ing the PEDOT:PSS for 10 minutes at 130 °C. The ac- 
tive layer consisting of poly (3-hexyl thiophcne) (P3HT) 
as donor and [6,6]-phenyl-C61 butyric acid methyl ester 
(PCBM) as acceptor was spin coated from a chloroben- 
zene solution (ratio 1:0.7). Finally, a second annealing 
step at 130 °C for 10 minutes and the evaporation of 
the metal Ca/Al cathode completed the solar cell assem- 
bly. The materials were purchased from H.C. Starck (PE- 
DOT:PSS), Rieke Metals (P3HT) and Solenne (PCBM). 
The current-voltage characteristics were recorded by a 
Kcithley 237 source measurement unit under an artificial 
AM1.5g sun spectrum with a spectral radiance of 100 
mW/cm , simulated by a 300W xenon lamp adjusted 
by a mismatch factor. Layer thicknesses and the ITO 
roughness were determined by a Veeco Dektak 150 pro- 
filometer. 



B. Numerical methods 

In order to simulate current -voltage characteristics, we 
use the elliptic differential equation system of Poisson 
equation, 



OF 



1 



(p(x) -n(x)), 



dx £o e 
two charge transport equations for electrons and holes 



(1) 



J„ = q(nfj, n F + D n dn/dx) 
J p = q{p^ p F - D n dp/dx), 



and the continuity equations in steady state, 
q(G(x)-R(x)) 
-q(G(x)-R(x)). 



dx 
dJp 
dx 



(2) 
(3) 



(4) 
(5) 



Discretized and combined with models for charge car- 
rier generation G(x) and recombination R(x), this can 
be solved numerically usi ng a pproximations according to 
Gummcl and ScharfetterPlHS] 

As basic parameters information on the charge carrier 
mobilities for electrons /i n and holes /x p , the correspond- 
ing effective charge carrier densities N c , N v , the temper- 
ature T and the material effective dielectric permeability 
e r are required. F denotes the electical field. 

Inside the blend, excitons generated by light on 
the polymer chains are splitted up on a femtosecond 
timescale when reaching a donor-acceptor (polymer- 
fullerene) interface.^ The hole remains on the donor's 
highest occupied molecular orbital (HOMO) whereas the 
electron is transferred to the lowest unoccupied molecu- 
lar orbital (LUMO) of the acceptor due to the energetic 
advantage of this state. Assuming an effective medium 



approach this process is allowed everywhere inside the ac- 
tive layer leading to free charge carriers which are trans- 
ported through the device. The energy difference be- 
tween the acceptor's LUMO and the donnor's HOMO 
energy is used as the effective band gap E g . For all 
presented calculations, the charge carrier generation rate 
G(x) = Gq is set to be constant throughout the device. 
The charge carrier recombination process inside the ac- 
tive layer which is bimolecular and non-geminate is de- 
scribed by the Langevin theory^l22l 



R = (/i„ + ix p ) (np - nf) 

eoe r 



(6) 



This rate is proportional to the product of charge carrier 
densities for electrons n and holes p until the thermally 
activated ground state is reached, which is indicated by 
the intrinsic charge carrier density 



m = y/N e N v exp (E g /2k B T) , 



(7) 



using the Boltzmann constant k B and the dielectric vac- 
uum permeability eo. 

According to the thermionic emission theory, we cal- 
culate the thermally activated charge carrier densities at 
the contacts by the energetic distance between the anode 
work function and hole conducting HOMO level at the 
interface, defined as the injection barrier <& p 



nth = N c exp {-E g + %)/k B T) 
p th = N v exp {-%/k B T) . 



(8) 
(9) 



Vice versa the injection barrier $„ is the difference be- 
tween the electron conducting LUMO and the cathode 
Fermi level at the cathode. 

In analogy, the local charge carrier density inside the 
active layer can be calculated from the energetic distance 
of conductive bands and their corresponding quasi- Fermi 
levels Ep n , E Fp , 

n = N c &tp((E Fn - E LUMO )/k B T) (10) 

p = N v exp((E H OMO-E Fp )/k B T). (11) 

In order to represent metal-semiconductor charge 
transfer processes in a general way, we use the surface 
recombination rate for electrons and holes defined by 



J s = qS(n~n th ). 



(12) 



Since this equation does not exclusively describe recom- 
bination, but an extraction of charge carriers out of the 
bulk it also can be considered as an extraction rate. In 
contrast to the thermionic emission theory the amount 
of surface charges n is changed by the surface recombi- 
nation velocity S, which defines the current transferred 
accross the metal-semiconductor interface. As we con- 
sider an active material between two metallic electrodes, 
we assume four independent surface recombination ve- 
locities S®'p, for electrons and holes at each electrode, as 
listed in Table [TJ Due to the different amount of charge 



carriers, holes and electrons at the anode are denoted 
as majority and minority charge carrieres. Vice versa at 
the cathode, holes and electrons are the minorities and 
majorities. Straightforward we distinguish between ma- 
jority and minority surface recombination velocities. 

All presented calculations were performed with the pa- 
rameter set given in Table [I] unless otherwise noted. 



parameter value 



description 



Ey 




1.1 eV 


$n, 


*p 


0.0 eV, 0.1 eV 


fJtn, 


Up 


10~ 8 m 2 /Vs 


L 




100 nm 


Go 




6.0- 10 27 m-V 1 


T 




300 K 


N c , 


N v 


10 26 m" 3 


e r 




3.4 



effective band garPEU 

injection barriers 

mobilities^ 

active layer thickness 

generation rate 

temperature 

effective density of states 

relative static permittivity 3 -^ 



surface recombination velocities: 



Op 

s: 



Op 



10 



CO 
-9 , 



i/s 



minority, cathode 
minority, anode 
majority, cathode 
majority, anode 



TABLE I. Standard parameters used in the macroscopic sim- 
ulation except as noted otherwise. An infinite surface recom- 
bination velocity is represented as 10 50 m/s in the numeric 
calculations. 
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FIG. 1. (Color online) Measured current-voltage characteris- 
tics of a bulk heterojunction solar cells without PEDOT:PSS 
layer under illumination. The cell with an additional 10 sec- 
ond oxygen plasma treatment of the ITO substrate shows a 
regular solar cell behavior whereas an extension leads to a 
distinct s-shaped behavior reducing the solar cell efficiency. 



prove this value. Longer treatment durations result in 
an equivalent s-shaped characteristics behavior but with 
a constant V oc level (not shown). 



III. RESULTS 



B. Numerical results 



A. Experimental results 

In Figure [l] typical current -voltage characteristics of 
the P3HT:PCBM organic solar cells without a PE- 
DOT:PSS layer under various oxygen plasma treatment 
timescales are shown. Within the first 10 seconds the 
power conversion efficiency increases from 0.72 % to 
2.9 %. The open circuit voltage is enhanced from 270 mV 
to 560 mV, the fill factor from 0.36 to 0.56 and the short 
circuit current from 7.37 mA/cm 2 to 9.23 mA/cm 2 . 

After 500 seconds of oxygen plasma treated the sam- 
ples undergo an efficiency loss to 0.0017 %. This can be 
attribute to the drastically drop of short circuit current 
to 0.047 mA/cm . The open circuit voltage decreases 
slightly to 235 mV, the fill factor to 0.15. The meassured 
roughness of the ITO substrate was not changed during 
the oxygen plasma etching within the profilometer noise 
level of 5 nm. The active layer thickness of all measured 
samples was about 230 nm. 

Deposing a PEDOT:PSS layer on top of the treated 
ITO anode slightly improves the initial solar cell work- 
ing efficiency. In contrast to a cell without PEDOT:PSS 
the first ten seconds of oxygen plasma etch do not im- 



In order to explain the experiment qualitatively, we 
focus on the effect of metal-semiconductor interfaces, 
namely a reduction of the surface recombination veloci- 
ties at the hole extracting contact (anode) for majority 
charge carriers. By limiting the majority surface veloc- 
ity from infinity (numerically rj 10 50 m/s) to 10~ 4 m/s, 
10~ 6 m/s, 10~ 9 m/s and 10 -11 m/s, s-shaped IV-curves 
are produced as shown in Figure [2j With the parameters 
given in Tab. IT] and infinite surface recombination veloc- 
ities, we retrieve an open circuit voltage of 609 mV. For 
the given reduced surface recombination velocities, V oc 
drops down to 360 mV. The calculated short circuit cur- 
rent decreases from 7.41 mA/cm to 4.34 mA/cm , the 
fill factor from 0.41 to 0.18. 



IV. DISCUSSION 

Metal-organic junctions posses an intrinsic charge car- 
rier density described by the injection barrier according 
to the thermionic emission theory. By this parameter no 
prediction on the amount of conducted charge is made. If 
there are any restrictions on the charge extraction process 
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FIG. 2. (Color online) Simulated s-shaped current-voltage 
characteristics due to a reduced majority surface recombina- 
tion velocity 5™ for holes at the hole conducting anode. With 
decreasing velocity, the s-shape gets more distinctive limit- 
ing the power converion efficiency by diminishing open circuit 
voltage, fill factor and short circuit current. At the lower right 
corner a magnification of the V oc dependence on decreasing 
surface recombination velocity is shown. 



this needs to be expressed by the surface recombination 
velocities according to Eqn. ( 12 1. If we assume a constant 



current that has to be transported through the interface, 
a reduction of the surface recombination velocity S has 
to be compensated by additional charges at the junction 
(n — n t h)i maintaining the passed currrent. These piled- 
up charges therefore create an applied voltage dependent 
space charge region which modifies the device energetic 
band structure. In Figure [2] the effect of constant major- 
ity surface recombination velocities at the anode is shown 
creating a typical s-shaped current -voltage curve. 

From literature we known planar semiconductor het- 
erojunctions to cause equivalent s-shaped characteristics. 
The physical origin of this effect is depicted in the upper 
row of Figure[3j The heterojunction is defined (I a) by the 
energy differences between the two LUMO ($«,) and the 
two HOMO bands ($ p ) of both semiconductors directly 
at the interface. Driven by an electric field (I b) addi- 
tional electrons above the intrinsic neutrality level can 
cross such an interface without transport problem if they 
are transfered from the energetically unfavorable level to 
a more favorable one. However, by band bending of the 
two semiconductors a local potential well is created if the 
polarity of the applied voltage is changed (I c) , charging 
a depending local space charg e and creating an s-shaped 
current-voltage char act erist ic EHE2 The same effect has 
been seen in copper indium gallium diselenide (CIGS) de- 
vices where such a s-shape behavior was assigned to by 
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FIG. 3. (Color online) Space charge accumulation in semi- 
conductor devices. (Upper row, a) At a semiconductor het- 
erojunction the two LUMO energy levels possess a work func- 
tion difference of $ n at the interface, (b) Additional electrons 
above the intrinsic concentration, injected by an external elec- 
tric field cross this interface from an electric unfavorable state 
to the more favorable one without restrictions, (c) If electrons 
are forced to overcome this energetic barrier in the opposite 
direction, a local potential well is established creating a pile 
up of electrons at the heterojunction. (Lower row, a) At a 
metal-organic interface the work function difference between 
the metal Fermi level and the semiconductor LUMO level $„ 
defines the intrinsic electron concentration at the interface. 
Due to the constant metal work function no potential well or 
space charge can be created, (b) By a finite surface recombi- 
nation velocity S n electron transport through the interface is 
limited. Electrons which are transported faster towards the 
interface than they are removed create a space charge, (c) 
The amount of injected charge carriers is reduced by finite 
surface recombination velocities creating a local charge car- 
rier depletion zone at the interface. 



light filled traps in an intermediate cadmium sulfide layer 
creating an illumination dependent energetic barrierP 

In device geometries such as bulk heterojunctions, 
without any planar semiconductor heterojunction, this 
approach cannot be applied. Due to the constant metal- 
lic work function metals are not able to create potential 
wells as shown in Figure^ (II). The injection barrier $„ 
(II a) defines the intrinsic charge carrier density at the 
interface as defined by the Boltzmann equation (Eqn. M) . 
Hence creating a local space charge needs the definition 
of an additional parameter, the surface recombination ve- 
locity. If, under extraction conditions (II b) charges are 
transported faster towards the interface than they can 
be extracted, those charges will pile up creating a space 



charge depending upon the applied electric field. Under 
injection conditions (II c) the injection of charges gets 
reduced if the surface recombination velocity is expected 
to be constant. Hereby the influence of the intrinsic holes 
rises, creating a local depletion zone. The injection bar- 
rier at metal-organic interfaces is not able to create s- 
shapes on its own. 

As explained in Figure [3] for electrons the issue holds 
for holes as well. For moderatly reduced values of 5° 
(10~ 4 m/s) (Figure^]), an s-shaped current-voltage char- 
acteristic is created at higher injecting voltages. It does 
not influence the solar cell efficiency parameters V oc , J sc 
and the fill factor significantly. Decreasing the surface 
velocity further, the deformation will approach V oc cre- 
ating a horizontal current plateau. Hence, all mentioned 
solar cell parameters are reduced. 



A. Energetic structure 

The deformation of the current -voltage characteristics 
originates from the internal energetic band structure and 
charge carrier distributions of electrons and holes. Fig- 
ureHlshows the device band structure under applied volt- 
ages of (a) V a = 300 mV and (b) V a = 900 mV, hence 
under extracting and injecting conditions, respectively. 
The energy levels of LUMO and HOMO bands (E LUM o, 
Ehomo) are indicated by the thick solid lines, those of 
quasi- Fermi levels E Fril E Fp by the blue and red dashed 
lines. At the anode (left contact) a limiting surface re- 
combination velocity of S p = 10 -9 m/s is assumed. For 
comparison, the underlying gray lines show the same so- 
lar cell with an infinite surface recombination velocity 
which will not lead to an s-shapc deformation. 

From these diagrams, position dependent information 
on the electric field, the charge carrier density and the 
current flux inside the device can be determined. The 
electric field is proportional to the slope of the energetic 
bands E LUMO and Ehomo {F = —qdE/dx). The lo- 
cal charge carrier densities are exponential functions of 
the energetic distance between the band energy level and 
the corresponding quasi-Fermi level (Eqns. ( fl0| ), (111). 
Multiplied with the slope of the related quasi-Fermi level 
these charge carrier densities result in the local current 
densities (J n = \i. n ndEp n jdx) caused by drift as well as 
diffusionPS 

Within the solar cell working regime charges gener- 
ated by light are extracted from the active layer into the 
metal electrodes. A reduced majority surface recombi- 
nation velocity S p will hinder an efficient holes extrac- 
tion by the anode and holes are piled up at the interface 
(Figure [41(a)). Consequently, at the anode interface the 
quasi-Fermi level of holes is shifted towards its HOMO 
level creating a quasi-Fermi level discontinuity between 
Ep p and the anode Fermi level. The interface hole den- 
sity is higher than predicted by the thermionic emission 
theory. Additionally to the already high hole concentra- 
tion at the anode interface for infinite surface recombina- 
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FIG. 4. (Color online) Energy structure of a solar cell with 
the energetic levels of HOMO and LUMO (solid and dotted 
lines) as well as their corresponding quasi-Fermi levels Ep n 
and Ef p (dashed red and blue lines) under illumination. The 
thick lines indicate the structure under a surface recombina- 
tion velocity of Sp — 10~ 9 m/s. Thin lines correspond to 
perfectly conducing contacts (Sp — > oo). In the case of hole 
injection by the anode (compare to Fig. wl (II b)) at an ap- 
plied voltage of 300 mV the quasi-Fermi level of holes split 
up at the anode. Describing the charge carrier density this 
generates a space charge combined with a higher electric field 
(oc dELUMo/dx). At a voltage of 900 mV holes are injected 
into the bulk by the anode (compare to Fig. [3] (II c)). Cre- 
ated by the surface recombination, less majority charges are 
present at the interface creating an electron depletion zone 
which decreases the local electric field. The cathode energy 
level (right side) was set to a fixed value of -4.25 eV. 



tion velocities additional majority charges will be added 
creating a space charge. Due to Coulombic interaction 
these charges increase the local electric field at the con- 
tact and consequently reduce the potential drop inside 
the semiconductor bulk. 3 -S In a solar cell dominated by 
drift, charges are transported more slowly through the 
device, thus increasing the probability of recombination. 
At voltages above the built-in potential additional 
charges will be injected into the semiconductor bulk 
by the metal contacts (Figure Kb)). A finite majority 
surface recombination velocity therefore will reduce the 
number of injected holes at the anode interface. Thus, 
the quasi-Fermi level for holes Ep p at the interface is 
shifted away from Ehomo- Within the interface region 
fewer holes and more electrons (due to the limited bulk 
recombination rate, Eqn. ffih) will be available. Conse- 
quently the electric field at the interface is reduced and — 
as in case of Figure U\h) — is even able to change signs. 
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FIG. 5. (Color online) Device band diagram at the open cir- 
cuit voltage for surface recombination rates of (a) infinity and 
(b) Sp = 10~ 9 m/s. Due to the condition of zero net current 
(J = 3n + jp — 0) holes pile up at the anode creating a space 
charge which directly reduces V oc - Whereas the quasi- Fermi 
distribution changes only slightly due to bulk recombination 
effects, but splits up directly at the contact the band bending 
and hence the electric field at the anode increases. 



The transition between both cases creates an s-shaped 
kink in the current -voltage characteristics of organic so- 
lar cells and leads to a space charge limited current at 
higher voltages. 



B. Open circuit voltage 

With decreasing surface recombination velocity the 
open circuit voltage drops according to Figure [2] In Fig- 
ure p^ the band diagrams for (a) infinite and (b) finite 
(Sp = 10~ 9 m/s) surface recombination velocities with 
their corresponding open circuit voltages of 0.61 V and 
0.47 V are shown. As in the case of charge carrier ex- 
traction the electric field at the anode is enhanced. Nev- 
ertheless the quasi-Fermi level distributions change only 
slightly except for their split-up directly at the anode 
contact. 

Under open circuit conditions, the total current at ev- 
ery location inside the device must be zero. Accordingly, 
the internal current densities for electrons and holes are 
equal and of different sign, J(x) — j n (x) + j P (x) = 0. 
Since an unlimited (infinite) minority surface recombi- 
nation always creates a steady diffusive recombination 
current towards the anode, this current has to be neu- 
tralized by an equivalent majority surface recombination 
current. If the surface recombination velocity according 
to Eq. ( 12 ) is reduced, the condition of zero net cur- 



rent can only be fullfilled by an accumulation of majority 
charges at the surface. The resulting creation of a space 
charge leads to an s-shape around V oc , and can even re- 
duce V oc (Figure [5]). 
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FIG. 6. (Color online) Open circuit potential plotted versus 
the hole surface recombination velocity at the anode. Assum- 
ing balanced charge carrier mobilities (fi n — fi p ) (a) shows 
a variation of the anode injection barrier from 0.0 to 0.4 eV 
in steps of O.leV whereas (b) presents the same solar cell un- 
der constant injection barriers and charge carrier mobilities of 
10" 10 m 2 /Vs, 10" 8 m 2 /Vs and 10~ 6 m 2 /Vs. At high surface 
recombination velocities no dependence of V oc from the sur- 
face recombination velocity can be observed in both cases. A 
reduction below a critical point results in a linear decreasing 
of V oc on the semi logarithmic scale which can be attributed 
to a logarithmic dependence of the local electric field at the 
contact on the majority charge carrier density. 



C. Analytical approximation of V oc 

Considering the dependence of V oc on the anode ma- 
jority recombination velocity (Figure p^ a) two different 
regimes can be distinguished. At high majority surface 
recombination velocities V oc remains constant whereas 
at lower values V oc decreases logarithmically. Addition- 
aly the open circuit voltage is lowered by the injection 
barrier $ p . 

The open circuit voltage can be calculated, using the 
local electric field F{x) — d/dx (qE H oMo(x)) bj^S 



t/ _ ^a ^ n ^p 

v oc 

q 



F(x)dx. 



(13) 



At high surface recombination velocities, V oc is domi- 
nated by the injection barriers. The barrier $ p changes 
the thermal equillibrium charge carrier densities n^, p t h 
(Eqns. ^, ([9])) to balanced conditions at the contact. 
Since the electric field generated by charge carriers is 
smaller for balanced densities, the band bending at the 
interface is reduced. Potential losses due to injection bar- 
riers are compensated, hardly reducing the effective open 



circuit voltage. If the band bending at the anode is fully 
compensated, any higher injection barrier therefore lin- 
early reduces V oc . 

With decreasing 5", additional majority charges are 
piled up at the anode. This enhances the electric field 
at the interface. In order to be able to influence V oc , the 
additional space charge must be significant compared to 
the electric field of infinite surface recombination solar 
cell contacts, a minimal and therefore critical additional 
charge carrier density is required, indicated by a kink on 
the V oc (Sp) diagram. 

In order to understand the slope of V oc on lower recom- 
bination velocities, the electric field at the contact can be 
calculated, solving the Poisson equation in one dimension 
(Eqn. (IT])) together with the current transport equations 
for electrons (Eqn. (pi)) and holes (Eqn. fcfy) for the elec- 
tric field. Under the condition of zero net current at V OCl 
we retrieve 



.,5 _ 



F{x) 



Dpdp/dx — D n dn/dx 
PH P + nji n 



(14) 



with the Einstein relations D n = finksT/q, D p = 
fipkgT/q for electrons and holes. As next to the anode 
contact 

H P p ^> fi n n and — fi p dp/dx ^> [i n dn/dx (15) 



are valid assumptions, the electric field simplifies to 



F = 



k R T 



(I 



Up) 



(16) 



in analogy to organic bilayer solar cells W^ Since the con- 
tact charge carrier density is the sum of thermally ac- 
tivated and accumulated charges by reduced surface re- 
combination {p — pth)i this equation can be transformed 
to 



q 



ln(p th ) + In 1 + 



P - Pth 
Pth 



(17) 



The left summand is the electrical field of a sample with 
infinite surface recombination. By the right side the ex- 
plicit effect of accumulated charges due to the majority 
surface recombination at the anode is considered. There- 



fore Eqn. ( 13 ) can be extended using the electrical field 
of a solar cell without blocking contacts F( s -kx>){z) and 
the separate charge carrier accumulation by the major- 
ity surface recombination velocity at the anode contact. 
Since at the cathode no additional holes are accumulated, 
we retrieve 



V or = 



E — <£> 



*, 



F t , 



k R T 



In 1 



)(x)dx 



■Is 



qSgpth 



(18) 



with a constant majority surface recombination velocity 
S®, a required surface recombination current J s by mi- 
nority charges and the thermally activated holes pth at 
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FIG. 7. (Color online) Calculated current-voltage character- 
istics with Voc set to voltage axis origin. The anode major- 
ity surface recombination velocity is reduced from infinity to 
10 -6 m/s. At voltages near V oc an ohmic transport regime 
(slope 1) can be found, which changes to a space charge lim- 
ited current at higher voltages (slope 2). The transition be- 
tween both transport regimes is indicated by an s-shaped de- 
formation of the calculated characteristics. 



the anode. An also reduced minority surface recombina- 
tion can limit J s and therefore prevent a decrease of V oc . 
By adding a PEDOT:PSS layer on the ITO substrate this 
effect was also seen experimentally. 

Beside of surface parameters, V oc is influenced by bulk 
recombinatio n pro cesses (Eqn. J6[)), hence the charge car- 
rier mobility! 38 * 39 ! In Figure ib) the V oc dependences 
on balanced charge carrier mobilities at 1CP 10 m 2 /Vs, 
10~ 8 m 2 /Vs and 10~ 6 m 2 Vs at fixed injection barriers 
of 0.1 eV and 0.3 eV are shown. High charge carrier 
mobilities possess a higher bulk recombination rate and 
therefore reduce V nr . 



D. Space charge limited current 

In order to extract information on the type of charge 
transport it is appropriate to set V oc as origin of the volt- 
age axis. Subtracting V oc as seen in the lower right sec- 
tion of Figure[2]from the voltage axis and choosing a dou- 
ble logarithmic presentation for the calculated current- 
voltage curves (Figure [7]) reveals an ohmic (J oc V) be- 
havior at low voltages. At higher voltages it asymptot- 
ically changes towards a trap free space charge limited 
current (SCLC) defined as 



9 V 2 

J SCLC = oM e oeryy- 



(19) 



This kind of charge transport was intensively studied in 
the pastp3E2 although the origin of s-shaped current- 
voltage kinks was not described. 

An s-shaped deformation created by a finite major- 
ity surface recombination velocity connects this behavior 
to the analytic SCLC prediction. Hence, in context of 
Fig. [71 an s-shape can be interpreted as a direct transition 
from an ohmic toward a space charge limited current. 



V. CONCLUSION 

Using an extended oxygen plasma etching process of 
the transparent ITO anode during organic solar cell 
processing, we were able to generate s-shaped current- 
voltage characteristics. Due to the absence of planar 
semiconductor hetero junctions we were able to exclude 
local potential wells as origin of the observed kink. Nev- 
ertheless, assuming finite surface recombination velocities 
this behavior was qualitatively reproduced in a numeric 
macroscopic device simulation. 

By analyzing the calculated energetic band structure, 
we assign space charges created by the reduced majority 
surface recombination velocity as origin of the character- 
istic deformation. Since also injection barriers change the 



charge carrier density at the contacts, we could see both 
parameters affecting V oc , but only the surface recombi- 
nation velocity being able to create an s-shaped current- 
voltage characteristic. Also in multilayer devices local 
space charges are responsible for the s-shape generation. 
We presented an analytic approximation showing that 
the loss of open circuit voltage on decreased majority sur- 
face recombination velocities can be countervailed by the 
reduction of the minority surface recombination veloc- 
ity. Finally we were able to show the s-shape indicating 
the transit from an ohmic conductivity to a space charge 
limited current. 
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